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SUMMARY 

Bounds for the tails of Dirichlet integrals are established by 
showing that each integral as a function of the limits is a Schur 
i unction. In particular, it is shown how these bounds apply to the 
simultaneous analysis of variance test and to the multinomial distri- 
bution. 
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MONOTONICITY PROPERT IES OF DIRICHLET INTEGRALS WITH APPLICATIONS 
TO THE MULTINOMIAL DISTRIBUTION AND THE ANOVA TEST 1 

Ingram Olkin 

Stanford University and Educational Testing Service 

1 . Introduction . The present paper is concerned with establishing 
bounds for the Dirichlet integrals: 

L l(a) = jf 1 ... / k fe t”" 1 (l-Zt.) d Ildt. , (1.1) 

0 0 1 1 

a a, , 

* 1 i* K lr vr 1 H ^ 

L 2 (a) = f ... J nt^ 1 (l-2t i ) d ndt. , (1.2) 
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(1.3) 



where a = (a,,..., a. ), and d is such that the integrals exist. 

The main result is that under certain conditions, -L^Ca), -LgCa), 
and -U(a) are Schur functions in (a. , • . • >a^.), (see e.g., Berge (1963) 
or Marshall, Olkin and Proschan (1967)). As a consequence of this fact 
it follows that if 



(a ^, « » • > 



a^) (b.^, • • • 
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in the sense that, after possible reordering, with a > ... > a , 

b, ^ ... ^ b, , that 
1 - - k 



m m k k 

Za • ^ Zb . , m =1, . . . ,k-l , Za » = Zb , 

I 1 ! 1 1 1 l 1 



(l.A) 



holds, then 



I^/a) < L^b) , L 2 (a) < L 2 (b) , U (a) < U(b) . 



(1.5) 



In this way we can generate many inequalities. Both the multinomial 
distribution and the simultaneous analysis of variance test are related 
to I> 2 (a) and U(a), respectively, and we provide several new results 
as a consequence of (1.4). 



2. The main results / We next prove the main results concerning 
L.. (a), Lg(a) and U(a). Theorem 1 is stated more generally and 
encompasses both L^(a) and L 2 (a). 

To show that a function F(a) s F(a 1 ,...,a k ) is a Schur function. 



we must show that 



, vV > 0 , 



( 2 . 1 ) 



for all i and j 



O 

ERIC 



4 . 



* 






Theorem 1. If f(x) is a non-negative monotone decreasing function, 



h rk k w -1 

(a;w) = / ... J f(2t) Ht. dt. 

J 0 Q 1 



( 2 . 2 ) 



where &1 > ... > a k > 0, 0 < W]L < ... < w R , then -L(a;w) is a 

Schur function in (a., , . . . ,a^) , whenever the integral is finite. 

Proof , Because of symmetry, we need only consider (i,j) * (l>2) in 
showing that (2.1) holds for -L(a;w). It is immediate from (2.2) -- 
noting the condition on f(x) -- that 



fJj(a;w) 



r2 r k k , V 1 k w -1 

■/ •••/ f <V|V a i fi dt i 



(2.3) 
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Condition (2.l) for -L(a;w) will be satisfied if 



/ ' 
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1 i k 

y’[a 1 x f( aj^-^gZ+Za t y. ) z 



w 2 ~l 



k w -1 

a 2 lf ^ a l Z+a 2 +2a i y i^ Z 1 ^ 
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(2.4) 



dz H dy <0 
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A sufficient condition for (2.4) to hold is that, pointwise, 



w -1 



a 2 f(a 1 -»a 2 z+Q)z 2 < a.^ (a^+a^z 1 



w, -1 



(2.5) 
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For 



O < z < 1, (aj+ajZ-K!) > (a^+a^! > s0 that f ( a l* a 2 z<< ^ - 
f (a^ta^) ■ Since 0.< V <y Z W 2 < z" 1 , and (2.4) holds. Il 

The result for U(a) is more delicate and a pointwise argument 
does not carry through. 



Theorem 2. 



k >°, 0 
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nt . 1 dt. 
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( 2 . 6 ) 



where d > Zw^ then -U(a;w) is a Schur function. 

Proof . Because of symmetry, we need only consider (i,j) = (l,2) xn 
showing that (2.1) holds for -U(ajw). It is immediate that 

du(a; 



w -1 k w.-l 
» a l Pi dt ^ 



a;w) = _ f f j- 

l l (l +a 1 +t 2 + **• +t k^ 



(2.7) 



Let t 2 -a 2 = z 
becomes 



and interchange order of integration; then (2*7) 



du( 

da. 



“ k w.-l 



and,, f ... f (nt , 1 dt,) f — L - 
1 { K 3 J 0 < 1+a l 



W -1 w -1 

« a. (z+a 2 ) dz 
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( 2 . 8 ) 
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We 



now use a pointwise argument on the inner integral with It^ 



fixed. Let z = (l +a 1 +a 2 +Zt i )v 5 SV then th ® inner integral becomes 



w -1 w^-1 



i 2 

“ a^ x (sv+a 2 ) 

J 0 s d_1 (l+v) d 



dv 



. ( du(a;w) d'j(a;w) < Q 

Consequently, n a ! - a 2> \ ^ ~ ca 2 / - 



provided 



w -1 

» (sv+a 0 ) 

^0 



Wg-i w 2 -l 
2 ' “ a 2 



V 1 



(sv+a. ) ] 

— dv > 0 . 



(1+v) 



The ordering 0 < w^^ < w 0 guarantees that the integrand be non- 
negative, so that (2.10) holds. || 



5, An application to the multinomial distribution . Let X = (X^ , 
have the multinomial distribution 



P{X=x) = ( n ) n ©j 1 

V' ,,X k 1 1 



where x = (x^...,**), 2x. =n, 0^. > ... > 6* > 0, 26. =1, and 

consider the tail probability P{X X > r,...,X fc > r l© 1 , • • • ,6 k ) , with 
r < n/k. 

Alam (1970) obtains lower and upper bounds for P{X 1 > r,...,X k 

9 ,...,©.} by averaging some of the 0’s, namely 
1 k 



(2.9) 



( 2 . 10 ) 
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(3.1) 



P{X 1 >r,...,X k >r!e*,e k ,...,O k ) < P{X 1 >r,...,X k >rle i ,...,0 k } (3.2) 

< P{X 1 >r > ... >X k >r 1 0. ♦...>? } , 

where. 0* = l-(k-l)0. and 0 = 10. A* It has been shown by Olkin and 

lc 1 

Sobel (1965) that P{X 1 > r , . . . > r 1©^ . . . ,0 k ) has a representation 

in terms of the Dirichlet integral 

/ 

I k 

where d«n-kr>0, and s(n,r,k) » r(n+l)/t( r(r)) r(n-kr+l)]. We may 

now make use of Theorem 1 with w^ = • • • *w k = r . Thus > if 0 < $ 1# 

0 .< q. < 1> and (p^,..»>P k ) ^ (<1^> • • • ><1^) > then 



>r,...,X k >rl0 1 ,...,0 k ) = K(n,r,k) / 1 



dt. 

l 



(3 



P{X 1 >r,...,X k >r|p 1 ,...,P k ) < P(X 1 >r,... / X k >rl 



I<1^> • • • ><1^) 



(3.4) 



The results of Alam are special cases of (3**0 since 




. >0fc) y (0^> • • • >0 k ) ^ (®> ,, *>®) 



Clearly, many other intermediate bounds can now be obtained. 

For the lower tail of the multinomial distribution (3*l)> we 
have the representation 
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C (0;r) s ?[x i < r,...,X m < i |0 1 ,..., e k ) A ( n , r , m ) (3.5) 



m m m-1 

l-ze, i-ze i -t 1 l- z t. 



- r 1 f yl — / 1 v«/frx 



e, e. 



m 



where m < min(k-l,n) and d=n-kr >0. As in Theorem 1, a direct 

differentiation of (3*5) with respect tc 0^, followed by the change of 

variables t -0 =v, t.-0.=z , yields 

2 2 J J J 



m-1 

e o «o’ v ®o' T ‘ l z i 



^-rr •••/ 

1 0 0 c 
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.{0. (v+e 9 )n(z +e ) f" 1 (e 0 -v-zz ) a drndz. ', 
x ^3 J j v j u ^ j 



m , m 
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where 0 Q = l-ZO^ By symmetry, 

m-1 

0. 0 -v 0 -v- Z z. 

ac ( 0 ; r ) r r r° 3 J 



d0. 
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m 



m 



0 0 



( )ii ( y « j > ) r ' 1 (» 0 - v -?3 )d av n dz 



3 - 5 m 



That -C:(0;r) is a Schur function follows from the fact that 

( e i" e 2^ ^ 8 i v+0 i e 2^ r 1 - -° for a11 v * 

As a consequence, we have the 



Corollary . If (p^,. ..,P fc ) (q^ • • *,< 1 ^), then 

P (X^ < r , . • • ,X m < r I , . . . , q^ ) > P (X^ < r , . . . ,X^ ^ r lp^ , . . . ,p^ ) 
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4 , An application to the simultaneous analysis of variance model . Suppose 
two hypotheses are tested using the same error variance for each test, so 
that we have 

4 A ■ V"! 

where q.^ and q 2 are independently distributed as X variates 

with n Q , n x , and n 2 d.f. respectively. Kimball (1951 ) obtained the 

inequality 

p{F l - F ia' F 2 - F zx ] - p(F l - F la ] P(F 2 - F aa ] ’ 



Where F and F are the 100 a percent points of the distributions 

la 2a 

of F^ and Fg. This inequality is of interest in that it provides a 
bound for the probability of making no errors of the first kind. We 
may use Theorem 2 to obtain a bound for P(F^ > F^, F 2 > Fgy) * 

Suppose that n^=n 2 =n, then 

» so ^n-1 jn-1 

dx dy , (4 , 2 ) 



p < p l- F a’ F fi - V * k / / 



C C (l+x+y) 



^n 0 -*n 



where c = nF a /n Q and k is a normalizing constant. Since (2c, 0) X(c>c), 
we obtain from Theorem 2 that 



P( F 1 > F a , F 2 > F a ) > P(F 1 > 2F a , F 2 > 0} = P{F 1 > 2F a ) . (4 . 3) 
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Since 



p{J iSV F a- V • P( W *W*SV ,p l f i^--'a ! - 1 ’ <4 ' 4) 



we obtain an alternative inequality to that of (4.l)> namely, 



P(F 1 < F a> F 2 < F a ) > 2P(F 1 < ¥ q ) - P(F 1 < 2F a ) 



(4.5) 



This is to be compared with 



F 2 <y >(P(F 1 <F 0 ))‘ 



( 4 . 6 ) 



We wish the larged bound, so that we need to determine the sign of 

[P(F 1 <F a )] 2 -2P(F 1 <F 0 )tl-P(F 1 >2F 0 )-[PIF 1 >F a )) 2 - p (F ;L >2F 0 ) . <4 - 7) 

It turns out that the difference is not always of one sign. When n=2, 

-n 

[P{F > c)] 2 = (1 +~) ° < (l+j^) 2 = P(F > 2c) , 

n 0 0 

so that (4.5) yields a better bound than (4,6) • With n Q •* », the 
difference becomes 



[P{X 2 > C )J 2 " p ( x n - 2cJ * 



(4.8) 
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m-1 . 

When n=2m, P{X 2 > c) = 2 c 3 /*'., and a straightforward analysis 

n 0 

shows that (4.8) is nonnegative. 

2 

For small values of c (less than c^n^ng))* [P{F>c}] is 

larger than P{F > 2c), whereas for c > c Q (n,n 0 ), P{F > 2c) is 
larger than [P{F > c) ) 2 , where c 0 (n,n 0 ) depends on n and n Q . 

As either n or n Q increase, the constant c 0 (n,n Q ) tends to 
increase. Since c will, in general, be of moderate size, it appears 
that (4.5) is the better bound in practise. 
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